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Lanczos Methods for UV- Suppressed Fermions 
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In this talk I indroduce lattice fermions with suppressed cutoff modes. Then I present Lanczos based methods 
for stochastic evaluation of the fermion determinant. 



1. INTRODUCTION 

There has been recent interest in the so-called 
'Ultraviolet slowing down' of fermionic simula- 
tions in lattice QCD [1-5]. These studies try 
to address algortmically large fluctuations of the 
high end modes of the fermion determinant. The 
goal is to increase the signal-to-noise ratio of the 
infrared modes and to accelerate fermion simula- 
tions as well. 

In fact, all the computational effort needed to 
treat UV-modes by above algorithms can be re- 
duced to zero by suppressing them in the first 
place [6]. The lattice Dirac operator of this 
fermion theory is given by: 

"I. ;,D\, , , , ^ 



D = -T 5 tanh ■ 



5^W/S/0 



/' 



where D w / s / is the input lattice Dirac operator, 
a the lattice spacing and /i > is a dimension- 
less parameter. For Wilson (W) and overlap (o) 
fermions as the input theory one has T$ = 75 . For 
staggered fermions is a diagonal matrix with 
entries +1/ — 1 on even/odd lattice sites. The 
theory converges to the input theory in the con- 
timuum limit and is local and unitary as shown 
in detail in [6] . The input theory is also recovered 
in the limit /i — > 00. For fi — > one has D — > /i, 
i.e. a quenched theory. 

Perturbative calculations with this theory are 
straightforward. In the following Wilson fermions 
are used as input. The inverse fermion propagator 
is given by: 



D(p) = —75 tanh — — 

a fi 



(2) 



with p = \j) Vl v = 1,... ,4} being the four- 
momentum vector. As usual, gauge fields 
are parametrized by su(3) elements, U(x) u = 
e iagA(x) v wit jj A{x) v g su(3) algebra, and the 
Wilson operator is written as a sum of the free 



and interaction terms: D 



+ D Z W . The 



splitting of the lattice Dirac operator is written 
in the same form: 

D = D° + D I , £° = ^ 75 tanh^ 
a fi 



where the interaction term has to be determined. 
This can be done by expanding D in terms of a//x: 

^ ^ r-n ,aHw 9 aHu/ .a 

d = d w [i + ci( — -f + C2 ( — -r + ■■■ 



(3) 

where C\,C2, ... are real expansion coefficients. 
Calculation of D 1 is an easy task if one stays with 
a finite number of terms in the right hand side of 
(3). Also, the number of terms can be minimized 
using a Chebyshev approximation for the hyper- 
bolic tangent. 1 

2. COMPUTATIONAL METHODS 

The effective action of the theory defined above 
can be written as: 



S cS = trf(A) 



(4) 



where A £ (£ NxN anc j j( s ) j s a rea i anc j smooth 
function of s G R + . The matrix A is assumed 



J I would like to thank Joachim Hcin for discussions related 
to lattice perturbation theory. 
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to be Hermitian and positive definite. Since the 
trace is dificult to obtain one can use the stochas- 
tic method of [7] . The method is based on evalu- 
ations of many bilinear forms of the type: 



T{b,A) = b T f{A)b 



(5) 



where b £ R w is a random vector. The trace 
is estimated as an average over many bilinear 
forms. A confidence interval can be computed 
as described in detail in [7]. Here I will describe 
an alternative derivation of the Lanczos algorithm 
[8], which is used in [7] and in [9] as well. The 
full details of this study can be found at [10]. It 
uses familiar tools (in lattice simulations) such as 
sparse matrix invertions and Pade approximats. 

The Pade approximation of a smooth and 
bounded function /(.) in an interval can be ex- 
pressed as a partial fraction expansion: 



/(»)«£■ 



fe=l 



dk 



(6) 



with Cfc £ WL, dk > 0, k = 1, . . . , to. It is assumed 
that the right hand side converges to the left hand 
side when the number of partial fractions becomes 
large enough. For the bilinear form I obtain: 



Ck 



fe=i 



A + d k t 



(7) 



A first algorithm may be derived at this point. 
Having computed the partial fraction coefficients 
one can use a multi-shift iterative solver of [11] 
to evaluate the right hand side (7). The problem 
with this method is that one needs to store a large 
number of vectors that is proportional to m. This 
could be prohibitive if to is say larger than 10. 
However, it is easy to save memory if only the 
bilinear form (5) is needed (see [10] for details). 

If a Pade approximation is not sufficient or dif- 
ficult to obtain, the Lanczos method is the only 
alternative to evaluate exactly the bilinear forms 
of type (5). To see how this is realized I assume 
that the linear system Ax = b is solved to the de- 
sired accuracy using the Lanczos algorithm. The 
algorithm produces the coefficients of the Lanczos 
matrix T n , which is symmetric and tridiagonal. 



Its eigenvalues, the so called Ritz values, tend to 
approximate the extreme eigenvalues of the orig- 
inal matrix A. In the application considered here 
one can show that [10]: 



E 6 ' 
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A + d k l 
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k=l 



ei 
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where 1„ £ M nx ™ is the identity matrix and e\ is 
its first column. From this result and the conver- 
gence of the partial fractions to the matrix func- 
tion /(.), it is clear that: 



F(b, A)^T n (b, A) 



|2„T 



ei f(T n ) ei (9) 



Note that the evaluation of the right hand side 
is a much easier task than the evaluation of the 
right hand side of (5) . A straightforward method 
is the spectral decomposition of the symmetric 
and tridiagonal matrix T n — Z n Q n Z^ , where 



ft n £ 



is a diagonal matrix of eigenvalues 



u)\, . . . , U! n of T n and Z n £ R" xn is the corrspond- 
ing matrix of eigenvectors, i.e. Z n = [zi, . . . ,z n ]. 
From (9) and the spectral form T n — Z n n n Z^ it 
is easy to show that: 



|2„T 



ei Z n f(n n )Z^ ei 



(10) 



where the function /(.) is now evaluated at indi- 
vidual eigenvalues of the tridiagonal matrix T n . 
The eigenvalues and eigenvectors of a symmet- 
ric and tridiagonal matrix can be computed by 
the QL method with implicit shifts [12]. The 
method has an 0(n 3 ) complexity. Fortunately, 
one can compute (10) with only an 0(n 2 ) com- 
plexity. Closer inspection of eq. (10) shows that 
besides the eigenvalues, only the first elements of 
the eigenvectors are needed: 



fn(b, A) 



l 2 £4/(^) 



(11) 



It is easy to see that the QL method delivers the 
eigenvalues and first elements of the eigenvectors 
with 0(n 2 ) complexity. 2 A similar formula (11) 
is suggested by [7] based on quadrature rules and 
Lanczos polynomials. Clearly, the Lanczos al- 

2 1 thank Alan Irving for the related comment on the QL 
implementation in [12]. 



Algorithm 1 Lanczos algorithm to compute (5). 

Set p = 0, pi = I/H&H2, qo = 0, qi = pib 
for 1 = 1,... do 

v = Aqi 
cti = qjv 

v := v - q l a i - qi-iPi-i 
Pi = \\vh 

qi+i = v/Pi 

Pi+l = -{PiOii + Pi-lPi-i) I 'fa 
if l/\pi+i\ < e then 

n = i 

stop 
end if 
end for 

Set (T n )i^i = Oii, (Tn)i+l,i = {T n )-Li+1 = Pi, 

otherwise (T n )ij = 

Compute uoi and zu by the QL method 

Evaluate (5) using (11) 



gorithm, Algorithm 1 has 0(nN) + 0(n 2 ) com- 
plexity, but it delivers an exact evaluation of (5) 
and for typical applications in lattice QCD the 
0(n/N) overhead is small. The method of [7] 
computes the relative differences of (11) between 
two successive Lanczos steps and stops if they 
don't decrease below a given accuracy. In order 
to perform the test their algorithm needs to com- 
pute the eigenvalues of Ti at each Lanczos step 
i. This is illustrated in Fig. 1 where the conver- 
gence of the bilinear form (5) and extreme eigen- 
values of A are plotted. The figure sugests that 
one can compute the bilinear form less frequently 
than proposed by [7]. 
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